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Abstract 

We consider the detection problem of a two-dimensional function from noisy 
observations of its integrals over lines. We study both rate and sharp asymptotics 
for the error probabilities in the minimax setup. By construction, the derived tests 
are non-adaptive. We also construct a minimax rate-optimal adaptive test of rather 
simple structure. 
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1 Introduction 

The problem of tomography is to reconstruct a two-dimensional function (image) from 
its Radon transform, i.e., from observations of its integrals over lines. This problem, 
and its extension to higher dimensions, appears in different scientific fields such as radio 
astronomy and medical imaging (see, e.g., [7J, [9], [20]). We consider the tomography 
problem from a statistical perspective that can be formulated as a problem of recon- 
structing a two-dimensional function from its noisy Radon transform (see, e.g., [5], [6], 
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Despite some work on the minimax estimation problem of a two-dimensional func- 
tion from its noisy Radon transform (see [6], [16], p2]), to the best of our knowledge, 
there exist no work on the corresponding minimax detection problem. The general 
statement of this problem is given in Section 2, while some preliminaries and notation 
in the minimax signal detection framework are presented in Section 3. Within this 
framework, in Section 4, we consider the detection problem of a two-dimensional func- 
tion from its noisy Radon transform and study both rate and sharp asymptotics for the 
error probabilities. By construction, the derived tests are non-adaptive. A rate-optimal 
adaptive test of rather simple structure is also constructed. The proofs are given in the 
Appendix. 



2 Formulation of the problem 

2.1 The Radon transform 

Denote by || • || the standard Euclidean norm in M 2 , i.e., ||x|| = (x 2 + x 2 ,) 1 / 2 , x = 
(xi,X2) G K 2 . Let H = {x G R 2 : ||x|| < 1}, be the unit disk in M 2 , and let \i denote 
the Lebesgue measure in R 2 . Consider the integrals of a function / : H i— >■ M over 
all lines that intersect H. The lines are parameterized by the length u G [0, 1] of the 
perpendicular from the origin to the line and by the orientation ip G [0, 2tt) of this 
perpendicular. Suppose that / G L^iH, fj) n L 2 (H,fi). Define the Radon transform of 
the function / by 



1Zf(u, ip) = — ^=^= / f(ucosip — t sin tp,u sin tp + t cosy?) dt, (u, ip) G S, 

2V1 - u 2 J-VT^ 

(2.1) 

where 

S = {(u,p) :uG [0,1], ip G [0,2tt)}. 

Thus, the Radon transform TZf is tt times the average of / over the line segment 
(parametrized by (u, cp)) that intersects H. It is natural to consider TZf as an element 
of L 2 (S,fj,o), where yUo is the measure on S defined by 

dfj, (u, ip) = U dip, (u, p) G S. 

TT 



2.2 The Gaussian white noise model 

Consider now the Gaussian white noise model 



dY £ (u,ip) = TZf{u,ip)dudp + sdW(u,ip), (u,p) G S, (2.2) 

where TV is a standard Wiener sheet on S (i.e., the primitive of white noise on S) and 
e > is a small parameter (the noise level). Although this model is continuous and 
real data are typically discretely sampled, versions of it have been extensively studied 
in the nonparametric literature and are considered as idealized models that provide, 
subject to some limitations, approximations to many, sampled-data, nonparametric 
models (see, e.g., 0, 0|, 0, TO- 

The Gaussian white noise model (|2.2p may also seem initially rather remote. One 
may, however, be helped by the observation that what it really means is the following: 
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for any function g G L 2 (S,fio), the integral 



g(u,<p)Kf(u,(p) dudap 

can be observed with Gaussian error having zero mean and variance equal to 
e 2 ff s g 2 (u,ip)dudip (see, e.g., @]). 

The Radon transform 1Z is a compact operator and its singular value decomposition 
(SVD) is well-known (see, e.g., [20]). To introduce it, let N = {1, 2, . . .} be the set of 
natural numbers, set Z + = NU{0}, and define a set of double indices giving rise to the 
following lattice quadrant 

r = {v:v = (j,l), j,leZ+}. (2.3) 
An orthonormal complex- valued basis for L 2 (H,fi) is given by 

Mr,6) = tt-^U + / + l) 1 / 2 ^ 1 (r)exp«j " 1)0}, ^T, (2.4) 

where x = (r cos 0, r sin 0) G H, with Z h a denoting the Zernike polynomial of degree a 
and order b, with a, b G Z + (see, e.g., [?]). The corresponding orthonormal complex- 
valued basis in L 2 (S,fio) is 

-i/v(u, ip) = ir~ 1/2 U j+ i(u) exp{i(j - l)ip}, v G T, (u, <p) G 5, (2.5) 

where 

sin((m + 1)0) 

£/ m (cos 0) = . - , m G Z+, G 0, 2tt), 

sm 

are the Chebyshev polynomials of the second kind. We then have (see, e.g., [6]) 

1Z((> U = b u ^ u 

with singular values 

K = vr(j + / + l)- 1 / 2 , i/er. (2.6) 

Since we work with real- valued functions /, the complex- valued bases (12.41) and 
()2.5p are identified, in standard fashion, with the equivalent real-valued orthonormal 
bases <\> v and ip v , v G T, respectively, defined by 

V2Re(4> u ), ifj>Z, 

if J = Z, (2.7) 
V2Im(^), if i < Z, 

with an analogous expression for ip v , v G T. 

Hence, by standard calculations (see, e.g., [6], [13]) and an application of the spec- 
tral theorem for the self-adjoint compact operator 7Z*1Z (TZ* being the adjoint of TZ), 
the Gaussian white noise model (|2.2p generates the following equivalent discrete obser- 
vational model in the Fourier domain, called the Gaussian sequence model, 

Vu = b v v + eS v , veT, (2.8) 

where y v = (TZf, ij) v ), v G T, are the "observations", b u , v G T, are the singular values of 
the Radon operator TZ, given by (|2.6p . V = (f,<f>v), v G T, are the Fourier coefficients 
of / with respect to 4> u , given by (|2.7p . and £„, v G T, are independent and identically 

distributed (iid) standard Gaussian random variables, i.e., £ u *~ M(0, 1), v G T. 
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2.3 The class of functions 



Crucial to the suggested detection methodology is the idea of considering minimax 
detection over certain classes of functions in / € L 2 (H, //) . Following [6] , we consider 
a special class of functions with polynomially decreasing coefficients 9 = {0„}„ e r, i-e., 
for some p > 0, L > 0, 



with 



e(p,L) 



eel 2 -. Yl V + 1 ) 2p ( / + x ) 2p ^ ^ l2 



(2.9) 



(2.10) 



It has been shown that J-(p, L) can be identified with the set of functions / which 
have 2p weak derivatives (provided 2p is an integer) that are squared-integrable on H 
with respect to the modified dominating measure 

dfi 2 p+i{x) = (1 - |M| 2 ) 2p dfi{x), x e H. 

This is weaker than the square-integrability with respect to /U assumed for the usual 
Sobolev spaces (see Proposition 2.2 in [T4"]). 



2.4 The aim 

The goal is to determine whether the two-dimensional function / corresponds to a 
known "etalon" function /o (i.e., to test the null hypothesis Hq : f = fo) or there exists 
a difference between / and /o (i.e., against the alternative hypothesis Hi : / = /q + A/ 
with A/ G F(p,L) (see (I2.9p - (l2.10p ). based on the observation of a trajectory {Y £ = 
Y £ (u,(p)}, (u,ip) G S 1 , from the Guassian white noise model (|2.2p . 

From mathematical point of view, we can take /o = by passing to the observation 
Y £ with dy e = dy e (u,0) — IZfodudcf). For this reason, without loss of generality, we 
assume in the sequel that fo = 0, use / in place of A/, and take the observation Y £ . 
In order to avoid having a trivial power (see below), our ultimate goal is to determine 
whether / satisfies f)3. If) (see below), using only tests calibrated in such a way that if one 
had run them in the absence of an / e J~(p, L), a certain restriction of the significance 
level (error probability) is met. 

In the sequel, we elaborate on the set under the alternative hypothesis and the 
suggested test statistics that provide a good quality of testing in the minimax frame- 
work. Before going into the details, however, we give the necessary preliminaries on 
the minimax signal detection framework in the standard Gaussian white noise model 
which provide the avenue for developing the suggested detection methodology and de- 
riving theoretical results for detecting a two-dimensional function from its noisy Radon 
transform. 

Hereafter, the relation A £ ~ B £ means that A £ /B £ — > 1 as e — > while the relation 
A £ x B £ means that there exists absolute constants < c\ < C2 < oo and £o > small 
enough such that c\ < A £ j B £ < c 2 for < e < Eq- 
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3 Signal detection in the Gaussian sequence model: the 
minimax framework 



Consider the Gaussian sequence model (|2.8|) . In order to avoid having a trivial min- 
imax hypothesis testing problem (i.e., trivial power), one usually needs to remove a 
neighborhood around the functional parameter under the null hypothesis and to add 
some additional constraints, that are typically expressed in the form of some regularity 
conditions, such as constraints on the derivatives, of the unknown functional parameter 
of interest (see, e.g., [11], Sections 1.3-1.4). 

In view of the above observation, the main object of our study is the hypothesis 
testing problem 

H : 9 = versus H x : ^a*0j|<l, ^6l>r 2 , (3.1) 

where 9 = {9 u } ue r G I 2 -, a v > 0, v £ T, and r £ > 0, r £ — > 0, is a given family. It means 
that the set under the alternative corresponds to an ellipsoid of semi-axes l/a„, v ET, 
with an / 2 -ball of radius r £ removed. (Here, I 2 = : YlueT < °°} with T given by 

Consider now the sequence r\ = {r] u } v er with elements rj u = 6 v /a v where we set 
a v = l/b u , v G T. In view of (|2.6p . the sequence r] = {r/ u }uer G h and the Gaussian 
sequence model (I2.8|) takes the form 

yv = r}v + ££ v , v^T. (3.2) 

The hypothesis testing problem (|3.ip can now be written in the following equivalent 
form 

Hq : 7] = versus Hi : rj G B(r e ), (3.3) 
where the set under the alternative, i.e., Q(r E ), is determined by the constraints 

9 = {r? G I 2 : Y. a2 M ^ Q ^ = 9 : > r 2 }, (3.4) 

i.e., the set under the alternative corresponds to an ellipsoid of semi-axes l/{a v u v ), 
v G T, with an ellipsoid of semi-axes r E /a u , v G T, removed. 

We are therefore interesting in the minimax efficiency of the hypothesis testing 
problem ()3.3[) - f)3.4|) for a given family of sets £ = B(r e ) C I 2 . It is characterized 
by asymptotics, as e — > 0, of the minimax error probabilities in the problem at hand. 
Namely, for a (randomized) test ip (i.e., a measurable function of the observation y = 
{Vu}uer taking values in [0,1]), the null hypothesis is rejected with probability tp(y) 
and is accepted with probability 1 — ip(y). Let P e „ be the probability measure for the 
Gaussian sequence model (|3.2|) and denote by E £jrj the expectation over this probability 
measure. Let a E (ip) = E £> Qip be its type I error probability, and let P £ (@ £ ,ip) = 
sup^gg^ E £ ^{1 — ip) be its maximal type II error probability. We consider two criteria 
of asymptotic optimality: 

(1) The first one corresponds to the classical Neyman-Pearson criterion. For a G 
(0, 1), we set 

&(e e ,a)= inf &(e £ ,v). 

We call a family of tests ip £ , a asymptotically minimax if 

a £ (ip £ ,a) <a + o(l), /3 £ (0 £ ,V> £ ,a) = /3 £ {@ £ ,a) + o(l), 
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where o(l) is a family tending to zero; here, and in what follows, unless otherwise 
stated, all limits are taken as e — > 0. 

(2) The second one corresponds to the total error probabilities. Let "y £ (Q £ ,^) be 
the sum of the type I and the maximal type II error probabilities, and let j £ (Q e ) be 
the minimax total error probability, i.e., 

7£ (e £ ) = inf 7£ (e e ,vO, 

where the infimum is taken over all possible tests. We call a family of tests ip £ asymp- 
totically minimax if 

7£ (e £ ,v £ ) = 7 £ (e £ ) + o(i). 

It is known that (see, e.g., [11], Chapter 2) that 

P e (Q e ,a) € [0,1 -a], 7 £ (© £ )= inf (a + f3 e (Q e , a)) G [0, 1]. (3.5) 

a£(0,l) 

We consider the problems of rate and sharp asymptotics for the error probabilities 
in the minimax setup. The rate optimality problem corresponds to the study of the 
conditions for which 7 e (O e ) — > 1 and 7 £ (O e ) — > and, under the conditions of the last 
relation, to the construction of asymptotically minimax consistent families of tests ip £ , 
i.e., such that 7 e (0 e ,V'e) — > 0. 

We are interesting in a set © e of the form 

e = 0(r e ) ={,£9 : H > r e }, 

where C h is a given set, | • | is some norm in I2 (not necessarily the standard I2- 
norm) and r £ — > is a given positive-valued family. For this case, we use the notation 
7 e (0(r £ )) = 7 e (r e ), /3 £ (0(r £ ),a) = f3 £ (r £ ,a) and we are interesting in the minimal 
decreasing rates for the sequence r £ such that J £ (r £ ) — >• 0. Namely, we say that the 
positive sequence r* — > is a separation rate, if 

7 e (r £ ) -t 1, and /3 e (r e , a) -)■ 1 - a for any a 6 (0, 1), as r £ /r* -)• 0, (3.6) 

and 

7 e (r e ) — )• 0, and /3 e (r £ , a) — )• for any a E (0, 1), as r £ /r* — > 00. (3.7) 

In other words, it means that, for small e, one can detect all sequences r/ £ 0(r e ) if 
the ratio r e /r* is large, whereas if this ratio is small then it is impossible to distinguish 
between the null and the alternative hypothesis, with small minimax total error prob- 
ability. Hence, the rate optimality problem corresponds to finding the separation rates 
r* and to constructing asymptotically minimax consistent families of tests. 

On the other hand, the sharp optimality problem corresponds to the study of the 
asymptotics of the quantities /3 £ (0 e , a), 7 e (0 e ) (up to vanishing terms) and to the con- 
struction of asymptotically minimax families of tests ip £ , a and ip £ , respectively. Often, 
the sharp asymptotics are of Gaussian type, i.e., 

Pe(& £ ,a) = $(H^ -u £ ) + o(l), 7 £ (0 £ ) = 2*(-« e /2) + o(l), (3.8) 

where $ is the standard Gaussian distribution function, is its (1 — a)-quantile, 

i.e., &(H( a ^) = 1 — a. The quantity u £ = u £ (r £ ) is the value of the specific extreme 
problem (|4.ip on the sequence space l 2 , and the extreme sequence of this problem 
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determines the structure of the asymptotically minimax families of tests Ve,a an d V'e- 
Moreover, we shall see that if u £ {r £ ) — > oo, then 7 e (r e ) — > 0, f3 £ (r £ ,a) — > 0, and if 
u £ (r £ ) — > 0, then j £ (r £ ) — > 1, (3 £ (r £ ,a) — > 1 — a, for any a 6 (0,1), i.e., the family 
u £ {r £ ) characterizes distinguishability in the testing problem. The separation rates r* 
are usually determined by the relation u £ (r*) x 1 (see, e.g., [TO], [H]). Hence, sharp 
and rate optimality problems correspond to the study of the extreme problem (|4.ip and 
of the asymptotics of the family u £ (r £ ). 



4 Minimax image detection from noisy tomographic data 

4.1 A general result: rate and sharp asymptotics 



Recall the Gaussian sequence model (|3.2j) . We are interested in the hypothesis testing 
problem (|3.3p with the set under the alternative £ = 0(r £ ) given by (13. 4p . 
Consider now the extreme problem 



2e 4 uee(r«)^ eJ . 



Suppose that 0(r e ) / and u £ > 0, and let there exist an extreme sequence {fj u } u( zr in 
the extreme problem (|4.ip , (Observe the uniqueness of a nonnegative extreme sequence 
{t^}, i/ST, because, by passing to the sequence {z v } u€l y with elements z v = fjl, 
vfT, we obtain the minimization problem of a strictly convex function under linear 
constraints.) Denote 

w u = —^== ^ =f , ^f, w = supw„, (4.2) 
V 2 2^er ^ ^ r 

and consider the following families of test statistics and tests 

te = Yl W » ( - l ) > = % e >H}, (4-3) 

where Q{A} denotes the indicator function of a set A. (Note that the values of fj u , w u , 
v € r, and depend on e, i.e., fj u = fj ViE , w v = w Uj£ , v E V, and wq = w$. £ .) 

The key tool for the study of the above mentioned hypothesis testing problem is 
the following general theorem. Its proof follows along the lines of the proof of Theorem 
4.1 in 1131 ; hence, it is omitted. 



Theorem 4.1 Consider the Gaussian sequence model h3.2\j and the hypothesis testing 
problem \3. 3\) with the set under the alternative given by l[3.4\ ). Let u £ be determined 
by the extreme problem \4- lj )> let the coefficients w u , v 6 T, and wq be as in and 
consider the family tests ip £t H given by |^.3| ). Then 

(1) (a) If u £ — > 0, then f3 £ (r £ ,a) — >• 1 — a for any a £ (0,1) and 7 e (r £ ) — > 1, i.e., 
minimax testing is impossible. Ifu £ = 0(1), then liminf f3 £ (r £ , a) > for any a £ (0, 1) 
and liminf 7 e (r e ) > 0, i.e., minimax consistent testing is impossible. 

(b) If u £ x 1 and wq = o(l), then the family of tests ip £t H of the form |^.3| ) with 
H = and H = u £ /2 are asymptotically minimax, i.e., 

a e(4> £ ,H(<*)) ^ a + o(l), 

P £ (Q(r £ ),i^ £;H{a) ) = (3 £ (r £ ,a)+o(l), 
7e (0(r £ ),^, Ue / 2 ) = 7e (n) + o(l) ; 
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and the sharp asymptotics \3. 8\) hold true, i.e., 

Mrs, a) = $(# (Q) - u £ ) + o(l), 
7 £ (r £ ) = 2*(-« g /2) + o(l). 

(2) If u e —7- oo, i/ien the family of tests t^ £ ,H of the form with H = T £ are 

asymptotically minimax consistent for any c G (0,1) and a family T £ ~ cu £ , i.e., 
7£ (e(r e ),V £ ,T £ )^0. 

Theorem 14.11 shows that the asymptotics of the quality of testing is determined by 
the asymptotics of values u £ of the the extreme problem (14. lh . In order to make use 
of it, one needs to study the extreme problem (|4.ip . This problem is studied by using 
Lagrange multipliers. Then, the extreme sequence in the above mentioned extreme 
problem is of the form 

1)1 = z 2 a 2 u (l - Aal)+, v G T, (4.4) 

where (i)+ = max(t, 0), t G R, and the quantities zq = zq j£ and A = A £ are determined 
by the equations 



E„ G r a ln, 



2^2 _ 2 



v • £ ) 



The equations (14. 5p are immediately rewritten in the form 

\ l = zlA-\j 2 , 



(4.5) 



(4.6) 



and, hence, the extreme problem (|4.ip takes the form 

u 2 e = e- 4 zlJ /2, (4.7) 

where 

Jl = ^2at(l-Aa 2 u ) + , 

J 2 = a£<#t£(1-A£) + , 

Jo = Ji- J 2 = J>2(1 - Aa^.. 

It is also convenient to rewrite (I4.6P and (14. 7p in the form 

Remark 4.1 Let u £ = ii £ (r £ ) be the value of the extreme problem (|4.ip with sequences 
a = {ajyjjygr and a = {a u } ut =r associated with the set under the alternative £ = 0(r £ ) 
given by (|3.4|) . and let u £ = u £ (r £ ) be the corresponding value of the extreme problem 
similar to (|4.ip with sequences a = Ca = {Ca u } u£ r and a = Da = {Da u } ue r hi 
(|3.4p . for some constants C, D > 0. Then, it is easily seen that the relation u £ (r £ ) = 
(CD)- 2 u £ (Cr £ ) holds true. 
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Remark 4.2 In order to obtain the corresponding rate and sharp asymptotics for 
the noisy tomographic data, we need to study the asymptotics of the quantities Jj, 
i = 0,1,2, given above. We note, however, that the methods used in [13] to study 
analogous asymptotics in a wide-range of linear statistical ill-posed inverse problems 
cannot be adopted to the problem at hand. The reason is that there does not exist a 
common ordering for the sequences a = {a v } u ^r and a = {<J u }u&r associated with the 
set under the alternative Q £ = G(r £ ) given by f|3.4[) . The arguments and techniques 
used to prove Theorems 14.21 and 14.31 below are specifically developed to tackle this 
problem. 

4.2 Rate and sharp asymptotics for the noisy tomographic data 

According to (|2.6p and (|2.1U|) . consider the double- index sequences 

a v = L-\j + l)P{l + l)P, ueT, (4.9) 
<Jv = vr-^j + Z + l) 1 / 2 , i/eT, (4.10) 

for some p > with T given by ()2.3|) . 



Theorem 4.2 Consider the Gaussian sequence model \3.%fy and the hypothesis testing 
problem Ii3. 3]) with the set under the alternative given by ft3.4\ )- Let {a u } u ^r and {ovji/er 
be defined as in |^.5| ) and |^.i0| ), respectively. Then 

(a) The sharp asymptotics \3. 8\) hold with the value u £ of the extreme problem j[ ) 
determined by 

On-I-I / 1 \ 1+3 /( 2 p) 
ul ~ 7T 4 L-^r 4 £ +3/p e~ 4 ( — — ) , (4.11) 

where B = Y%Li m_3 > L2021 < B < 1-2022. 

(b) The asymptotically minimax family of tests i^ £ ,h <ire determined by the family 
of test statistics t £ given by \4.3\ ) with coefficients w u , fGT, and wq as in ^4.2\ ), and 

3 rl 



4p+3 



with extreme sequence {fj u }uer satisfying \4-$ with A 
(c) The separation rates are of the form 

r * =£ 4p/(4 P+ 3)_ (412) 

Remark 4.3 It is easily seen that the asymptotic results in Theorem 14.21 hold true 
uniformly over p £ £ for any compact set £ C (0, oo). 

Remark 4.4 Rate and sharp asymptotics in the corresponding minimax estimation 
problem under the L 2 -risk have been obtained in |6j. In particular, the asymptotical 
(as e — > 0) minimax rates of estimation are given by 

i? £ 2 :=inf sup E\\f-f\\^e^^\ 
f f£F(p,L) 

where the infimum is taken over all possible estimators / of / based on observations 
from the Gaussian white noise model (|2.2p . (Here, we adopt standard notation and 
write g\ (e) X 52 (e) to denote < liminf ((71(e) /52(e)) < lim sup (51(e) /52(e)) < 00 as 
e — > 0.) By comparing r* with R £ , it is observed that the asymptotical minimax rates of 
testing are faster than the corresponding asymptotical minimax rates of estimation; this 
phenomenon is common in nonparametric statistical inference (see, e.g., {llj . Sections 
2.10 and 3.5.1, 
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4.3 Adaptivity and rate optimality for the noisy tomographic data 

The family of tests considered in Section 14.21 depends on the parameter p that is usually 
unknown in practice. Therefore, it is of paramount importance to construct families of 
tests that do not depend on the unknown parameter p and, at the same time, provide 
the best possible asymptotical minimax efficiency. Such families are called adaptive (to 
the parameter p) and the formal setting is as follows. 

Let E be a compact set in (0,oo) and a family r £ (p), p G E, be given, where 
e > is small. Let the set @ £ (p,r £ (p)) be determined by the constraints (|3.4p with 
&v = o, v {p), v £ I\ and r £ = r £ (p), and set 

e £ (E) = \Jq £ {pMp))- 

pes 

We are interesting in the following hypothesis testing problem 

Hq : 7] = 0, versus H% : rj G £ (E). 

We aim at finding conditions for either 7 £ (0 £ (E)) — > 1 or 7 £ (6 £ (E)) — > 0, and to 
constructing asymptotically minimax adaptive consistent families of tests ip £ d such 
that 7 £ (6 e (S),^ d ) ^ as 7e(© e (S)) ->• 0. 

Let tf £ (p) = u £ {p,r £ {p)) be the value of the extreme problem (j4. 1 j) for the set 
Q £ = e £ (p,r £ (p)). Set 

u £ (E) = inf u £ (p). 

We are interesting in how large u £ (E) should be in order to provide the relation 
7 £ (@ £ (E)) — > 0. We say that the family u £ d = u £ d (E) — > oo characterizes adaptive 
distinguishability if there exist constants < d = d{Y>) < -D(E) = D < oo such that 

7e(©e( s )) -> 1 as limsup pgS n £ (p)/n^ < d, 
7£ (G £ (S))^0 as liminf peE n £ (p)/< d > 

Observe that it follows from the asymptotics (|4.1ip that, by making r £ (p) larger or 
smaller, one can increase or decrease u £ {p,r £ (p)) in order to get u £ (p,r £ (p)) ~ u £ , for 
all p 6 E and any family u £ > 0. 

We call a family r^ d (p), p £ E, such that x it £ (p, r £ d (p)), the family of adaptive 
separation rates. 

Note that the relation 7 £ (@ £ (E)) — > is possible if n £ (E) — > oo. However this 
implication does not hold for the tomography problem under consideration, as we show 
below. (A similar situation appears in some ill-posed inverse problems, see |13j).) 
Hence, hereafter, adaptive distinguishability conditions and adaptive separation rates 
are sought for the tomography problem. In contrast to Theorem 14.21 there is price to 
pay for the adaptation. We show below that 

uf = Vlogloge- 1 , (4.13) 

yielding a loss in the separation rates in terms of an extra factor ydog log e _1 in e. 
Furthermore, the derived families of tests are of simple structure. (A similar loss in the 
separation rates was first observed in [22] and more recently in [13j.) 

Specifically, let p be unknown, p G E, where E = \p m i n ,Pmax], < p m in < Pmax < 
oo, be a compact interval in (0,oo). Let us also consider the collections 

p fc GE, c k ~2/r £ (p k ), k = 0,l,...,K, K = K £ x logfe" 1 ) log log^ 1 ), 
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where p = p ma x > Pi > ■ ■ ■ > PK = Pmin (the collection p k , k = 1, 2, . . . , K - 1, will 
be specified in the proof) and collection of statistics i £iCfc of the form 

te,c k = Yl W ^{(yl/e 2 ) - 1), w v>k = ^ w 2 U)k = i (4.14) 

for v G T given by (|2.3[) and C v> k = {y • a UtPk < c&}. Consider the following families of 
thresholds and tests 



H £ = 20og(K £ ), y £ = {y : t e>Ck < H £ , V < k < K £ }, A = ^y e , (4.15) 

where A denotes the complement of a set A. 
Denote also 

Hp) = t~~t~o-> H^) = IHp) ■ v e e} c (o, oo). (4.16) 

4p + 3 

Theorem 4.3 Consider the Gaussian sequence model \3.2) and the hypothesis testing 
problem h3. 3]) with the set under the alternative given by \3.4ty - Let {a u } u ^r end {o- u } u ^r 
be defined as in |^.ff[ ) and |^.i0[ ) ; respectively. Then 

(a) (lower bounds) Let the set 0(E) given by j4-16\ ) contains an interval 
[a, b], < a < b < 4/3. There exists constant d = d(Ti) > such that if 
limsup pgE n e (p)/ v / loglog(e- 1 ) < d, then 7 e (0 E (E)) -)• I. 

(b) (upper bounds) For the family of tests ip £ given by ft4.15\ ), a{ip £ ) = o(l) and 
there exists constant D = -D(E) > such that if lim inf pe s u £ (p) / \J\og log(e^ 1 ) > D, 
tten/3 e (e 6 (E),^) = o(l). 

(c) (adaptive separation rates) The adaptive distinguishability family u £ d is given 
by !{4-L3\ ) and the adaptive separation rates r £ d (p), p 6 E, are given by 

4p/(4p+3) 



rf(p) = e^loglog( £ -i) 



Remark 4.5 Rate and sharp adaptation in the corresponding minimax estimation 
problem under the L 2 -risk have been obtained in [6]. In particular, [6] showed that 

Rl = p £ (p,L)(l + o(l)), 

where 

1 / 4 \ 2p/(2p+2) 

P ^ L) = 2 {wTzj) {{2p + 2 ) L ) 2/(2p+2) e4p/(2p+2) ' 

and constructed an adaptive penalized blockwise Stein- type estimator / (see [6], equa- 
tions (3.12) and (5.8)) such that 

y E\\f-f\\ 2 , 

lim sup sup p — — — = 1, 

for any < p\ < p2 < oo and < L\ < L2 < 00. Note that, in contrast to adaptive 
minimax separation rates, there is no price to pay for adaptation in the corresponding 
minimax estimation problem under the L 2 -risk (i.e., a global measure). The unavoid- 
able logarithmic factor for adaptivity that appears in the minimax separation rates 
r £ d {p) stated in Theorem 14.31 also appears in various other adaptivity problems, such 
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as minimax signal detection (see [13] ) and minimax estimation under the L 2 -risk (see 
[24] ) in some ill-posed inverse problems. It also resembles the minimal price one needs 
to pay for adaptation in minimax estimation under the Z 2 -risk (i.e., a local or pointwise 
measure) that has been observed in [18] , [2] , [19] and [23] in the case of Lipschitz and 
Sobolev balls and, more recently, in [3] and [T] in the case of Besov balls. 



5 Appendix: Proofs 

For simplicity in the calculations, we omit the factors L" 1 and tt~ 1 in (|4.9p and (I4.10p . 
respectively. In other words, from now onwards, we work with a u = (j + l) p (l + l) p 
and a v = (j + 1 + l) 1 / 2 , p > 0, v G T. The final results can be obtained on rescaling by 
using Remark 14.11 



5.1 Proof of Theorem 4-2 



It follows from Theorem 14.11 that the efficiency in the detection problem under consid- 
eration is determined by the asymptotics of the quantity 

2 fre\ A Jo 



where 



Ue ~ V E J 2 J 2 ' 



Ji = Ji(A) = ^^(l-^a 2 )+, 
J 2 = J 2 (A) = AY^al4(l-At£ 
Jo = J (A) = Y, a t(l-Aal) 2 + = J 1 -J 2 , 

for i/ £ V given by (]2.3p ). Moreover, the quantity A = A e — > is determined by the 
relation 

rl = A J -l. (5.1) 

In order to study the asymptotics of u e , we are interesting in the asymptotics of the 
functions Ji(A), i = 0, 1, 2, as A — )■ 0. We first, however, start with the asymptotics of 
the following function 

i(A) = Yl ^ veT - ( 5 - 2 ) 

{v. Aal<\} 

Proposition 5.1 Let 1(A) be defined as in $5.2\) . Then, as A — > 0, 

J( A )^^ A -3/(2 P ) 5 (53) 

3 

where B = m " 3 > 1 - 2020 < B < 1-2021. 

Proof. Set j + 1 = m, I + 1 = n, H = [m -1 ^ -1 /^] and H x = [A' 1 /^], where [t] is 
the integer part of t. Consider the set 

C m , n ,p,A = {(m, n) : m > 1, n > 1, (mn) 2p < A^ 1 }. 
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Then, we have 

1(A) = (m + n-l) 2 

(m,ra)eC m , niP ,A 

Hi H Hi Hi 

= 2 ^ ^2(m + n- l) 2 - ^ ^(m + n- 1) 

m=l n=l m=l rt=l 

Hi 



= - V ((m - 1 + H)(m + H)(2m - 1 + 2fT) - (m - l)m(2m - 1)) + 0(A -1/p ) 

m=l 
Hi 

= 2 ^ (Hm 2 + F 2 m - Hm - H 2 /2 + F 3 /3 + F/6) + 0(A~ 1 / p ). 

m=l 

We now indicate the asymptotics of the items in the last sum. Observe that 

H = -A- 1 /^ + a(A, m), H x = A' 1 '^ + /3(A), 
m 

where a(A,m) G [0, 1) and /3(A) € [0, 1). Thus, we have 

Hi Hi 

Y^Hm 2 = ^2(A- 1 /^m + a(A,m)m 2 ) = A- 1 / p /2 + 0(A^^), 

m=l m=l 

Y,H 2 m = f](A- 1 /(2p)/ m + a ( Am ))2 m ^^-i/P log(A -i/(4p) ) = A-VPlog^- 1 ) 

m=l m=l 
Hi Hi 

Y^Hm = Y( A ~ 1/(2P) +a(A,m)m) ~ A" 3 ^, 

m=l rra=l 

f> 2 /2 = ^(i-VW/^a^.jfxi- 1 * 

m=l m=l 
Hl 4~3/(2p) H i 

= ^-3— E-- 3 + 0(^), 

m=l rra=l 
Hi Hi 

= ^(i- 1/(2p) /m + a(Am))/6 = o(A- 1/p ). 

m=l m=l 



Therefore, 

1(A) = N m H + 1/p ) ~ — ^ 3/(2p) . 

m=l 

The proposition now follows. □ 

Let us now return to the asymptotics of Ji(A), i = 0, 1, 2, as A — > 0. Introduce the 
following function 

F(t)=I(r 1 )= Y (m + n-l) 2 , t>0, (5.4) 

(m,n)eC mj „ iJ(jt _i 
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and observe that F(t) is nondecreasing in t > such that F(0) = 0. It follows from 
Proposition 15.11 that 

F ( t )^^3/(2p) ) (55) 

For T = A -1 , the functions Ji(A), i = 0, 1, 2, could be rewritten in the form 



Ji(A) = / (l-t/T)dF(t), 
Jo 

MA) = f T (t/T-(t/T) 2 )dF(t), 
Jo 

Jo(A) = [ (l-t/T) 2 dF{t) = J l -J 2 . 
Jo 



Integrating now by parts, we get 

J X (A) = T- 1 [ F{t)dt. 
Jo 

In order to study the asymptotics of the above integral, we divide it into the following 
two parts 

r-T XT 1 / 2 pT 

/ F(t)dt = Si + S 2 , Si= F(t)dt, S 2 = F(t)dt. 

Jo Jo Jt 1 / 2 

Hence, it suffices to check that S\ = o{S 2 ) and to use the asymptotics (I5.5|) of the 
function F under the integral in S 2 . It is easily seen that 

< Si < F(T 1/2 )T 1/2 x r 3/(4p)+1/2 

and that 



Therefore, we get 



Wp - A-*/™. (5.6) 
u ' 3(2p + 3) v ' 



Similarly, for T = A , we get 

J 2 (A) = T~ l j tdF{t)-T- 2 j t 2 dF(t) =2T~ 2 J tF(t)dt 

J \3(4p + 3) 3(2p + 3) 



4BP -A~^ (5.7) 



(4p + 3)(2p + 3)' 



and that 



Thus, it follows from (|5TTj) . (i5T6|) . (loTTD and dHHD that 

A — r 2 U 2 ~ r 4+3/ % -4 2P±3 / 1+3/(2p) ( } 

4p + 3 e £ e 25 V 4 P + 3 
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Therefore, the separation rates r £ are of the form 

r * =£ 4 P /(4 P +3)_ (51Q) 

In order to get the sharp asymptotics, in view of Theorem 14.11 it is enough to check 
the condition 

max {l/:a 2 A<1} o-2(l _ Aa? v ) 
wq = —= 1 = = o(l), v G r. 

^E { „:aiA<i } vt(l-Aal) 2 
This condition follows directly from the relation 

r "K? 1 ^ = °( J 0( j4 ))- ( 5 - U ) 

{u:alA<l} 

Indeed, for m = j + 1 £ N, n = I + 1 g N, the last condition follows from (see 

at = (m + n - l) 2 < AA~ l ' p < A~ 3 / 2p x J (A) 

as A — )• 0, for v = (m — 1, n — 1) such that a 2 , = (mn) 2p < A^ 1 . 
The theorem now follows. 



5.2 Proo/ o/ Theorem \4^3 



Let p be unknown and consider p G E, where E = [fW„,;Wr], < p min < p maa; < oo, 
is a compact interval in (0,oo). Recall the expression <p{p) from (I4.16j) and that Z + = 
NU{0}. 

We first obtain the lower bounds. Take a collection p^, k G Z + , such that 

= a + k5 £ , k = 0, 1, . . . ,K = K E 

with 



K \og{e^Y 
Therefore, 

p fc G [b' 1 -3/4, a" 1 -3/4] = E, b < 4/3, 

and 

4 = <KpaO - HPk-l), Pk-l ~Pk = (4pfc + 3)(4p fc _i + 3)J e /16 X - — -3-— = o(l). 

log(£ !) 

Assume without loss of generality that, uniformly in p G X, 

u £ {p) ~ Vdloglog(e- 1 ), (5.12) 

where the constant <i > will be specified below. This corresponds to taking, uniformly 
in p G E, 

r e (p) ~ (e^loglog^ 1 )) 174 )^, (5.13) 

where 

2B \ -(2p+ 3 )/( 2 p) 



= da(p), a(p) = 2p + 3 \l<t>(p) 
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Take 

T fc ~ (2r e (p fc )r lM % k = 0,l,...,K e . 
By construction, we have 

T k -T k ^ ~ T fc _x (2^-i-^ 1 (e^dloglog(e- 1 ))-^ - 1 
= T fc _ 1 (exp(log(2)(l + o(l))) - 1) ~ T fc _!. 

Observe that the function F(i) = defined by (I5.4p depends on p. Also, by (|5.5I) . 

we can write 

F(t) = F p (t) = E (n + i 2 + l) 2 

{((n+i)(i 2 +i)) 2p <t, ix,i a e^+} 

= E ^ = ^i(t 1/p )~^ 3/(2p) . (5-14) 

{v. al<t) 

Set 

A fc = {^o = (J, G 1+ x Z+ : T fe _! < (j + l)(f + 1) < T fc }. 
Take a collection > 0, Ufc )1/0 , z^o G ^fcj = 1, 2, . . . , if, such that 



0, otherwise, 



with ^ UQ = ±1, vq E Ajfc. We then have 

E <*o< = 4 E < = «2 W^) - ml-x)) ~ = 2rl(p k ) (5.15) 

and 

E 4,^l <, Ph =4 E «, Pfc <2r fc 2pfc r £ 2 (p fc )(l + (l)) = i(l + (l))<l. 

(5.16) 

Furthermore, we have 

ul(pk) ~ d loglog(e _1 ), cZ > 0, p k G S. 
Consider now the priors 

1 K 

^fc = II (^fc^o e "o + 5 -zfe^o e -o )/ 2 ' fc = 1 ' 2 ' • • • ' K > n = ft E 7Ffe ' 
^oeAfc fc=i 

where {e l , } 1/[)g 2 +x z + is the standard basis in the space I 2 that corresponds to sequences 
indexed by £ *Z+ x Z+, and <^ is the Dirac mass at the point rj £ I 2 . The relations 
(|5TT5jl and (15TT61) imply 

7r fc (e £ (p fc ,r £ (p fc ))) = l, 7r(e e (E)) = l. 

Let P njs = E nk P e)V , P n = E^Psft be the mixtures over the priors. It suffices to check 
that 

E efl ((dPjdP^ - I) 2 ) = o(l). (5.17) 
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Using evaluations similar to Section 5.6 in [13], we have 

K 



E £fi {{dPjdP £fi - l) 2 ) = ^ J2 E £>0 (idP v JdP eja - if 

k=0 
1 K 

k=0 

^ h ( 2 £ S mh 2 (. fc 2 </2 £ 2 ) ) - 1 ) . 
fc=o \ \ i/oeA* / / 



5 A 

Note that <t 2 = j + I + 1 < Tf. if i^o G A&. Therefore, uniformly over uq £ A^, we have 

fM> < X ^2+2^-2 X e 2/(4 Pfc+ 3) (loglog(e -l )){ 2 Pfc+ 2)/(4 Pft+ 3) = q(1) 

and (since sinh 2 (zf cr^/^e 2 ) ~ 2|cr*,/4e 4 ) 

2 £ sinh 2 (z 2 a 2 /2 £ 2 ) ~ ^" 4 r e 2 (p fc ) ~ || 2^r 4 + 3 /^K 4 

~ ||2 3 /PMa(p fc )loglog( e - 1 ). 
One can take d > such that, for any p^ £ S, 

75 - 7 fe=o,i,...,K E 1 2p fc + 3 V 3 J J 

Then we have 

E £fi [(dPjdP £fi -l) 2 ) < ^Wexph £ sinh 2 (z 2 </2 £ 2 )] -l] 

fc=o \ \ i/oeAfc / / 

Klog^Hg' 1 ) log^He" 1 ) _ m 
if 2 log^" 1 ) ° UJ ' 

We now obtain the upper bounds. 

Similarly to the proof of the lower bounds, assume without loss of generality that 
u £ (p) ~ \J D log log(e _1 ), uniformly inp G S, where the constant D > will be specified 
below. This corresponds to (|5.13p with d replaced by D, uniformly in p £ S. 

In order to evaluate the type I error probability, we consider a different grid with 
different K = K £ , i.e., 

(p(p k ) = a + k5 £ , k = 0,l,...,K = K £ , <f>(p K ) = b > a > 0, 

where 

S = 8 E - (6 " a) log(2) 



K £ log(e x ) log log (e x )' 
Let us evaluate the exponential moments 

£? E) o(exp(/ii £)Cfc )), /i > 0. 



17 



Recall that (y v /e) *~ A/"(0, 1), v G T, under Po- Recall the set C^fc = : < c k }, 

f € r, A; = 0, 1, . . . , K. Let the family h = h £ be taken in such way that 

h max u^fc = o(l). 

Then, we have 



k 



= exp ^2 (~ hw u,k ~ log(l - 2hw v>k )/2) 

= exp ^2 h 2 wl tk (l + 0(hw U)k )) 

= exp(/i 2 /2)(l + o(l)) . (5.18) 
Let h = H £ . Then, for k = 0, 1, . . . , iT e , we have 



n max = /i £ -j- < 



' ' ' " V 2E^,X) 1/a (4< Pfc (l + o(l))/3) 1/2 



2 



x H £ r l J^\p k ) x e 2 /( 3 -Nf*)(] og i og ( e -l))2(l-H» fc )/(3-H Pfc ) = o(1 ) j 
and by (|5.18p . for any k = 0, 1, ... , _fQ, we have 

P £ ,o(k Cfc > ff.) < ^' 0( e e ^ ( ( fg) £ ' Cfc)) ~ exp(fi*/2 - H 2 ) = eM~H 2 £ /2) = K 
This implies that, for the type I error probability, 

aWe) < Y, P ^(ts,c k > H £ ) < K-\l + o(l)) -> 0. 

fc=0 

Let us evaluate the type II error probability for 

v g e e (s) = U e 6lP (r e (p)). 

pes 

There exists p such that r\ G £iP (r £ (p)), p k < p < p k -i- Observe that 

Pe(y,il>e) < min P £>ri (t £ ,c k < H e ). 
Denote /i £jCfe = E £t11 (t £ , Ck ). We then have 
^, Cft = e ~ 2 «fc,fe^; Var e>IJ (t eiC J = l + 4e" 2 ^ u; 2 , ^ = i + 0{h £)Ck ). (5.19) 

Let us now evaluate 
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where the function 

F P (c) = J2 °t 

{u:al p <c} 

is of the form (j5.14j) with asymptotics given by (|5.14p as well. Observe that a u , Pk = 
and, hence, 

E « = E « " E « * ^(p) - c- 2 ^ = r 2 (p) (l - — ') . 
veC„ >k v£F {v.a u , Pk >c k } \ r e\P) c k I 

Because a(p) p ^ p ^ 4: satisfy the Lipschitz continuity, we have by (|5.13p with d replaced 
by D, 

r £ {Pk) 
r £ (p) 

and 



(1 + 0(4)) exp((p0(p) - K </»(p fc ))(log(e- 1 ) + (log(Z)- 1 - logloglogfc- 1 ))/^), 



o < P m - Pk <t>{ Pk ) = -u(p k )-m) < -a ~ :{loft(2) 



4 vv ^" y vw^- 4 -6 4io g ( e -i)loglog(e- 1 )' 
We get 

1 < < 2 3 /( 41oslog ( £ " 1 »(l + o(l)) = 1 + o(l), (5.20) 

r e (p) 

and, moreover, we have 

P ^P a 1 

— = H , < p-pk = < p k _i -p k ' 



Pk Pk log(e 1 )loglog(e x ) 

Thus, we have 

c fc < c£ /pfc = c fc • ct p/pk =c k (l + o(l)), r £ (p)cf> fc = 2(1 + o(l)). 
These relations and ()5.20|) imply, for e > small enough, 

E > r 2 £ (Pk)(l - 1/4(1 + o(l)) = r 2 e (p k )(3/4 + o(l)). 

v&C v>k 

Therefore, for large enough D > 0, we have, for all k = 1, 2, ... , iC e , 

r £ 2 fa)(3/4 + o(l)) r £ 2 (p fc ) / 3^ A 

^ - JMM) ~^l^ +0(1) J 



( J Da(p fe )loglog( e - 1 )) 1 /2 2 -3/(2p fc ) fi^+ (l) 



> 2tf £ ~ 4Vloglog(F T ). (5.21) 

It follows from (|5.19p and (|5.2ip that one has, for k such that rj £ £jP (r £ (p)), p k < 
P < Pk-i, 



Pe,r}{te,c k < -^e) — Pe,r){h £ ,c k ~ te,c k > h £ ,c k ~ H £ ) < 



Var £jr? (/i £]Cfc ^£,c fc , 



< ^ " hl Ck < 41oglog( £ -i) +0? ' (1) - 0(1) - ^ 
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Hence, for any w £ 8 E (E), one has 

PeMs) < l/(41oglog( £ - 1 ))+ 0r; (l) =o(l). 
Pe(9e(2>),i>e) = SUp /3 £ (r]^ e ) < — 1 + SUp 0^(1), 

r,ee e (s) 41oglog(e L ) v ee e (S) 

where, by (I5T221) . (gJH) and (l5TT9j) for any n G £ (E) 

m i£ ~ 2 ^2ueC Uik w l,kVl 4max veCk w v , k 2 -1^-1/2 
°^ 1 ) = 7P2V VT^f- < — ~(loglog(e )) /. 

This completes the proof of part (b) of the theorem. 

Part (c) of the theorem follows immediately from parts (a) and (b) of the theorem 
and (j4.1ip . The theorem now follows. 
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